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Abstract 

We apply an asymmetric version of Kirman's herding model to volatile financial 
markets. In the relation between returns and agent concentration we use the square 
root law proposed by Zhang. This can be derived by extending the idea of a critical 
mean field theory suggested by Plerou et al. We show that this model is equivalent to 
the so called 3/2- model of stochastic volatility. The description of the unconditional 
distribution for the absolute returns is in good agreement with the DAX independent 
whether one uses the square root or a conventional linear relation. Only the statistic 
of extreme events prefers the former. The description of the autocorrelations are in 
much better agreement for the square root law. The volatility clusters are described 
by a scaling law for the distribution of returns conditional to the value at the 
previous day in good agreement with the data. 
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1 Introduction 



Any model for the price or index of a financial market should account for the 
following so called stylized facts [1,2]. (i): The sign of the return cannot be 
predicted. In most models this fact is build in by assuming the return pro- 
portional to an iid noise, (ii): The probability for having an absolute return 
larger than x decays for large x as a power law x~ nT . The tail index ny seems 
to be universal [3] in the range 3-4. A precise determination of rij- for dayly 
returns is hampered by the impossibility to reach the asymtotic regime in 
the data. Therefore one should require a good description of all returns, (iii): 
Returns are not independent. They form so called volatility cluster. A quan- 
titative measure of this effect may be the distribution of first passage times 
[4]. In this paper we use the distribution of the returns conditional to the 
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value of the previous day. (iv): Absolute returns are correlated. In the case of 
high frequency data a hyperbolic decay of the autocorrelation indicates a long 
memory process. For dayly returns it behaves neither as a power law nor a 
simple exponential expected from a Markov process. 

To acchieve these properties one can distinguish two schools. The first are the 
stochastic volatility models. Here one assumes the return to be an iid noise 
multiplied with a volatility factor \fVr-, where Vt typically follows a Markov 
process. Examples are the GARCH model [6] and its continuous time version 
[7]. Other posibilities are the CIR- model [8] or the 3/2 model [9] discussed in ref 
[10]. The multifractal model [11] implements fact (iv) by occurence of many 
time scales. Any model for dayly returns faces the following difficulty The 
quality of data for dayly returns does not allow to estimate reliably more than 
three parameters, for example a scale parameter given by the first moment, 
a combination of parameters describing the tail index and a rate constant for 
the time dependence as in GARCH(1,1). Therefore the parsimonity princi- 
ple invoked by Calvet and Fisher [11] is rather a must and not a freedom. 
For stochastic volatility models a second difficulty arrises, that the choice of 
the free parameters is dictated by mathematical convenience. Their relation 
to market properties remains obscure. This difficulty is avoided in agent or 
microscopic models (for a recent review see ref [12] and the references given 
therein). The behaviour of the agents may be characterized by utility functions 
[13]. The observed power law (fact (ii)) motivates the use of critical physical 
models as the application of percolation by Cont-Bouchaud [14]. The perco- 
lation cluster are interpreted as herding behaviour of the agents. In the Lux- 
Marchesi model [15] a herding mechanism appears explicitly in the transition 
probabilities between chartists and fundamentalistic traders. This leads to an 
intermittent behaviour of agent numbers and thereby also for the returns. To 
see whether herding alone can describe the stylized facts a symmetric Kirman 
model [16] has been studied [17]. Quantitative agreement with the data could 
be acchieved [18] within an asymmetric Kirman model. In the latter agents 
differ slightly from those in the Lux-Marchesi model that the two types of 
chartists (optimists and pessimists) are replaced by one sort of agents (noise 
traders). The difference between optimists and pessimists is described by an 
iid random number. With this approximation the return behaves similar to a 
stochastic volatility model, except \JV~t follows a Markov process instead of 
Vt- In this paper we apply the Kirman model with a further simplification. 
In empirical investigations it has been found [18,19] that fundamentalists are 
much less affected by herding than noise traders, which implicates a very 
asymmetric herding model. In the limit of large asymmetry the number of 
free parameters is reduced from three to two. In the spirit of parsimonity we 
adopt this approach. 

For an application to financial markets an equation of state is needed relating 
the return to ratios of agent numbers. In a Walrasian market [18] the volatil- 
ity factor \fVr depends linearely on the concentration of agents. The model 
using the simplified Kirman model and this linear relation has been studied in 



2 



[19]. In this paper we adopt another equation of state. In contrast to a Wal- 
rasian market the time for an agreement of the agents on a common return 
ought to be finite. Assuming a random walk Zhang [20,21] concluded that V t 
itself should be proportional to the agent concentration. We use this approach 
in the present paper. It offers two desirable properties. Since in stochastic 
volatility models Vt and in herding models the concentration proportional to 
Vt follow simple Markov processes, one might be able to show the equivalence 
of the latter with one of the stochastic volatility models. Secondly one can 
use the empirical observation of Plerou et al [22], that the return and the 
imbalanced volume satisfy the same relation as the order parameter and an 
external field in a critical mean field theory. In the framework of a universal 
$ 4 [23] there ought to be another equation for the return at zero field. We 
will show that Zhang's law emerges within such a model naturally from risk 
aversity. Whether the Walrasian linear relation or Zhang's law are correct can 
be decided at present only empirically. We will compare our model and the 
results from ref [19] with the probability density function (hereafter abbrevi- 
ated with pdf ) and the autocorrelation function with the empirical returns 
derived from the DAX. 

The paper is organized in the following way. In section 2 we describe our ver- 
sion of the Kirman model and in section 3 a derivation of Zhang's law. The 
resulting unconditional pdf and the common pdf for two absolute returns at 
different times are given in section 4. The square root law together with a 
Gaussian noise allows analytical calculation of various observables. As exam- 
ples we treat the autocorrelation function and the pdf conditional on the value 
of the previous day. These predictions we confront with the data in section 5. 
In section 6 we make some concluding remarks. 



2 Kirman's Herding Model 

We want to describe the changes of prices in a stock market with different 
agents which can change their strategy as function of time. Apart from the 
specification of an agent model a relation between the relative changes (return) 
of the stock index and the agent numbers is needed which will be treated in 
the next section. As herding model we use the Kirman model [16] discussed 
by Alfarano et al [18]. It is defined by the transition probability between two 
types of behaviour. S denotes the mood of noise traders and F a conservative 
or fundamentalistic strategy. The probabilities per unit time to change opinion 
are given by 

ir(S -> F) = (a 2 + b ■ n F ) 

7r(F^S) = (a 1 + b-n s ) . (1) 
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The parameters ai j2 correspond to spontaneous changes of mind. The second 
terms in equ (1) are proportional to the agent numbers ng,F in the new state 
and therefore describe the herding effect. Despite of the extensive character 
of 7r this Markov process leads for a large number of agents to a non trivial 
equilibrium distribution for intensive ratios as u oc rip / {ns + np). As shown in 
[18] the master equation belonging to (1) leads to a Fokker Planck equation 
(hereafter abreviated by FPE) for the pdf w(u,t) and a Langevin equation 
for u(t). In application of the model it turned out [18,19] that in most cases 
ai >> a 2 holds. This means that fundmentalists are much less influenced by 
herding than noise traders. Invoking the parsimonity principle we consider the 
limit of large a\jb. For the ratio 

a\ np 
u — 

o n F + n s 



(2) 



it is shown in appendix A that one obtains the following FPE 



dw(u,t) d 
dt du 



u 



d 

[3-1 + —U 
ou 



w(u, t) 



(3) 



with the parameter (5 = a 2 /6 — 1 describing the equilibrium pdf. This implies, 
that the number of fundamentalists is much smaller than that of the noise 
traders. If we transform the Langevin equation belonging to (3) to the variable 
V = l /u 

Ay = bAt y t (1 + (1 - 0)y t ) + (2bAt , (4) 



we obtain a well known model of stochastic volatility, the so called 3/2 model 
[9] provided we identify y t in the relation between the return r t and a Gaussian 
noise r\ t 

r t = yfVt ■ Vt (5) 



with the volatility factor V t . Equ (4) corresponds to a special version of the 
3/2 model where drift and diffusion constants are related to ensure existence 
of an equilibrium pdf for y. From V t = u t follows another frequently used 
model [8]. 

The limit of large a\/b affects also the time scales. As shown in the appendix 
the process (1) should be applied to micro time steps At m , whereas the variable 
u changes with the observed time scale At. Both are related by At = a\At m jb. 
Eliminating a±/b with (2) we get 



which shows an explicit dependency of the agent ratio on the time resolution 
At. 

As shown in [9] the FPE can be solved analytically for the equilibrium pdf 

w (u) 



w o( u ) = r ^ + ^ exp(-u) (7) 



and for the conditional pdf w(u, t + t \u, t ) = A(u\u, exp(-bt)) with 



, (uuz) 13/2 ( z(u + u)\ (2y/uuz\ 
A(u\u,z)=w (u) \ _ z exp f /^f ^L_j . (8) 



1/3 denotes the Bessel function with imaginary argument. The kernel A has 
the following convolution property 

J duA(u\u, z)A(u\u', z') = A(u\u', z ■ z') (9) 

and can be used for the time evolution of w(u,t + t ) 

w(u,t + t ) = J duA(u\u,exp(—bt))w(u,to) . (10) 



3 Derivation of Zhang's Law 



The problem of this section deals with relation of agent ratio to the observable 
return r t . The small decay time in the autocorrelation of the returns suggests 
that the dynamic of the price is much faster than the change of agent numbers. 
The extreme case of an instantaneous Walrasian market has been assumed in 
ref [18] to model the behaviour of agents in the Lux-Marchesi model[15]. With 
the additional assumption that the change between optimistic and pessimistic 
noise traders occurs also on a fast time scale one gets [18] a linear relation 
between the return and the agent ratio ns/np. Using u for small values of % 
one gets 

r t = A\n S t = r ■ — ■ r] t (11) 

The noise r] t incorporates many changes during At due to the noise traders. 
In this Walrasian approach the price change is set instantaneously, a closed 
order book is assumed and the agents have to accept any new return irre- 
spectively of their risk aversity. Considering for example the XETRA market 
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these assumptions seem to be questionable. Removing the first assumption 
Zhang [20] assumed the time needed to get the new price is proportional to 
the demand, which is linear in the agent numbers. Describing the evolution of 
the price by a random walk the return will be proportional to the square root 
of the agent numbers. Therefore the linear relation (11) should be replaced by 
the so called sqare root law given by 

r t = r - — U • rj t (12) 



It has been applied [21] in context of the percolation model of Cont-Bouchaud 
[14] in order to obtain a tail index larger than 2. Equ (12) agrees with the 
return obtained in the 3/2 stochastic model. Therefore the parameters of the 
latter can be related to the behaviour of agents. 

To obtain a less qualitative derivation of (12) and to include the effect of risk 
aversity and an unbalanced order book we start with an observation made by 
Plerou et al [22], who showed empirically that the absolute return v increases 
with the imbalanced volume Q at small time intervalls At and small Q as 

v t = \r t \ oc Q} t /S (13) 



with an exponent 5 = 2.9 ± 0.3. At large At and Q the return saturates as 
v = B tanh(Bi$Y). As suggested already in ref [22] this could be interpreted 
as an equation of state in a critical mean field theory if v corresponds to a 
scalar order parameter and Q to an external field. As dynamical fields in such 
a theory we use the return $ fc j expected by agent % of type k = F, S. After 
time At equilibrium is reached at a common return r = E^^]. A critical 
mean field theory with a scalar order parameter is described by the universal 
$ 4 model [23] . This is defined by the following cost or energy functional 

S(*) = S int (AQ) + E E 

k=S,F i 



k.i 



(14) 



The first term S[ n ^ is a positive definite quadratic form in the differences 
A^ikji = $fe,i — describing the interaction between the agents. We assume 
it strong enough that in equilibrium the mean field value for is independent 
of i, k. This means the agents reach a common value v for their expected 
absolute return. We normalize by 5^(0) = 0. The <3> 4 terms lead to large 
cost functions for large <£>. They account for the risk aversity of the agents. 
We expect that noise traders are much less affected by risk aversity than 
fundamentalists and can set /3s to zero. Since the disordered phase E[&] = 
is not observed, the coefficients in the quadratic term in (14) must be 
positive. A negative contribution to the cost function of the noise traders can 
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be interpreted as 'happy loser' effect [20], that fundamentalists are willing to 
accept a loss of money restricted by their risk aversity. The analogous term 
for the fundamentalists can be neglected by setting a F = 0. The last term 
in (14) corresponds to the coupling of the return to an imbalanced volume. 
If this coupling is independent of k the coefficient can be set without loss of 
generality to 1. In the limit of large agent numbers the model is solved by a 
constant <&k,i = v with negligeable variance. The value of v is given by the 
minimum of (14). With the above discussed approximations S reads 

S( v ) = ~r n F v 4 — 7^ n s v 2 — £lv(ns + n F ) ■ (15) 



Introducing the parameter Q by 



n = 



(nsasf 
(n s + n F ) 2 f3 F n F 



1/2 



(16) 



the value of v making S(v) to a minimum satifies the following mean field 
equation 

( n s®s V _ n s a s v ^L = q ( 17 ) 
\(n s + n F )£l J (n s + n F )Vl tt 



The solution of (17) can be written as 

v n s + n F 



fl a s n s V^o 



Equ (18) expresses the scaling property of the critical theory. Q, il and v 
may be rescaled by the same factor without altering the results. The scaling 
function / does not depend on any parameter. Its behaviour near x ~ and 
x — > oo is given by 

, 1 + \x + • • • for i«l 

f(x) = { 2 (19) 
(x) 1 / 3 (l + (x)- 2 ^ + ---) for x»l 



If only the leading terms are kept, for historical reasons the first form is called 
the zero field equation of state and the second the critical isotherm. The latter 
written in v reads as 

v3 = ^±^. tt (20) 

p F n F 
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where we recover (13) for Q 3> Oo- For Q <C Q we get another equation of 
state 



as n s , 
^ = A/75- • — - , 21 
V Pf n F 



which reproduces the square root law of Zhang. The size of Q depends on 
the time resolution At, if we use the Kirman model of the previous section. 
Inserting equ (6) for the agent ratios we get 

/An 1/2 

oc . (22) 



Therefore the empirical relation (13) should hold only at small time scales in 
agreement with the data [22]. For large time scales the Zhang's law ought to 
be applied. In addition we have the prediction, that the rate factor r in (12) 
should vary with At as 

r oc (At) 1/2 . (23) 



This power law is in nice agreement with the observed [30] power of 0.51. 
There are still two problems within our derivation of the square root law. The 
application of a critical theory is confined to small values of the fields. At finite 
values the scaling law (18) may no longer valid. Secondly the return does not 
saturate as required by the data, but will increase as (20) for large fl. There 
exist of course many models which have the same small field expansion as (14). 
Even if we restrict ourselves to a minimal model of Ising spins which has been 
used in ref [22] to describe the saturation there will be an extra parameter 
for the distance from the scaling law (18). As shown in appendix B for a spin 
model the mean field equation (17) should be replaced by 

ns^s v = 3 _ / n s a s vy _ ^ 

(n s + n F )Q (3-y 2 )y \(n s + n F )Q Q 3 tt J 



with a parameter y describing the deviation from the critical region. For y — > 
we recover the critical mean field equation. Finite y may be adjusted to the 
observed saturation of the return with the unbalanced volume Q. Therefore 
the critical isotherm (20) is very model dependent. On the other side the zero 
field equation for Q <C fl is almost unchanged. In the case of (24) we obtain 



v = aiy) . f.^ (25) 
V Pf n F 
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with a y-dependent constant a(y) which is the non zero solution of 
a-y= 3 ^_ y 2 ' tanh ( a • v) ■ 



(26) 



Since this constant can be absorbed into the rate constant r , the square root 
law appears to be robust against deviation from criticality. In appendix B we 
also replaced the Ising spin model by a compact model, which leads to the 
same conclusions. 

At present stage of research only empirical observations can decide on the 
possibilities (12) or (11). In the present paper we want to discuss whether 
there are advantages of the square root law (12) over the linear relation (11). 



4 Unconditional probability and 2-point function 



Suppose one observes a time serie of the absolute returns H t = [vi, ■ • ■ , v t ]. We 
assume the time t to be an integer by absorbing the time step At in the rate 
factor b. Only the combined process [H t , U t \ with the history U t = [u , ■ ■ • , u t ] 
of the herding variables is a Markov process. This property is expressed by 
the following recursion formula for the pdf f(H, U) 

f(H t , U t ) = g(v t \u t ) ■ A(u t \u t ^, e" 6 ) • f(H t _ u U t ^) . (27) 



Assuming the sqare root law (12) with a Gaussian noise the conditional dis- 
tribution g is given by 



v 2 u s 



^(u|«) = Aov^-exp(-- T ) (28) 







with the normalization constant N = ^2/irrQ. To calculate the pdf for the 
absolute returns H t the solution of the recursion (27) has to be integrated over 
all u variables. The sqare root law leads to an exponential dependence of g as 
function of u, which matches with the exponentials in A allowing analytical 
calculations. In principal the recursion (27) can be used to compute f(H t ). 
The singularity of A for exp(— b) ~ 1 presents a problem, if one encounters a 
small value of the rate b. For the determination of the parameters f3, r and b 
two simpler cases are sufficient and we defer an approximative solution of the 
recursion (27) to a future publication [24]. 

The common pdf f(v t+to ,v t ) for the absolute returns at t and t + t is ob- 
tained from the t times iterated recursion (27) by integrating over U t and all 
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intermediate Vt+i ■ ■ -Vt-i+t - The latter integration eliminates the correspond- 
ing g-factors and we can use the convolution property (9) of A to perform the 
intermediate u integration leading to 



f(v t+to ,v to ) = J du du to g(v t+to \u)A(u\u to ,e bt ) 

J dH^dU^fiH^Ut,) . (29) 

For large t the operator A converges to the equilibrium distribution wo(u) 
independent of u to . Therefore f(v t+to ,v to ) factorizes in the limit of large t as 

/ K+to , «to ) = G o (vt+t ) • / (v t0 ) (30) 



with the equilibrium pdf 

G (v) =Jdu g(v\u) w (u) = ^o^^ (l + ' ^ (31) 

This Pareto like distribution becomes at large v a power law v~ nT ~ 1 with a 
tail index nx determined by the herding parameter 

n T = 2(3 + 2 . (32) 



From equ (31) all equilibrium moments may be calculated ( see equ (D.2) in 
the appendix). For example the second moment 



(33) 



will be used for expressing r in terms of (3 and the observable second moment 
in equilibrium 1 . 

To derive an explicit form of the two point function f(v t+to ,v to ) we take in 
equ (29) large to with finite t. Analogue to the previous case the distribution 
factorizes and we obtain 

f(v t+to ,v to ) = J dudu'g(v t+ t \u)A(u\u,e~ bt )g(vt \u)wo(u') . (34) 

The integrations in (34) can be carried out and f(v t +t ,v to ) is expressed in 
terms of hypergeometric functions. The exact form is given in appendix C. An 
interesting property of / is that it depends only on a combination £(i> t+to , v to ) 

1 Expectation values taken with the equilibrium pdf (31) will be denoted by Eq[ . ] 
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but not each of the variables v to and v t +t separately. Such a behaviour implies 
scaling laws. We demonstrate this scaling law in the simple case of t = 1 
and neglected terms of order b. As shown in appendix C the conditional pdf 
f(v t+ i\v t ) is given by 



r(/3 + 2) , vf +1 



-13-2 



If we eliminate r with equ (33) and introduce the scaling variable x t +i by 

Vt+l 

/2(3E [v 2 ]+v. 



V t+1 / Q „N 

x t+1 = == , (36) 



we get 

f , | x 2 r(/? + 2) / 2 r /3-2 



The conditional distribution of x t+ i is independent of the previous return v t . 
From this we can determine /3, since only at the true f3 the data for x t+ i with 
different v t will collapse into a single scaling function f(x t +i). This function 
can be compared with the predicted curve (37). 

Knowing the two point function the autocorrelation of v q can be derived (see 
appendix D ). The ratio 

R = ^El^ = F{q/2 > q/2] P + 1; 6 ~ bt) (38) 



involves again a hypergeometric function. For negative values ofo; = g — (3 — 1 
the variance of v q exists and we find the autocorrelation function C(t) 

C<l{t) = vSR ' ( F(g/2 ' q/2] e ~ M) " X ) • (39) 



If we expand for large t the hypergeometric function we find with the leading 
term an exponentially decay 



For small b and finite t the deviation from the exponential dependence (40) 
and the exact expression (39) may be substantial. It is interesting to note that 
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Table 1 

(5 values and their statistical errors obtained by a fit in four different time intervalls. 
The last row gives the fit using all data. Column 3 gives the x 2 value per data point. 



for a > the ratio R becomes singular at t — 0. Applying Eulers relation to 
the hypergeometric function one observes for the ratio 



r(/? + i)r(q) 

R ~ iw ' m 



(41) 



a power law in t. This does not mean a long memory since neither the variance 
exists nor it consists in a large t effect. The example given by Lillo et al.[28] 
shows that Markov processes may lead to power law in t also with existing 
variance. 



5 Empirical Comparison with the DAX 



In this section we want to estimate the parameters of the model described 
in the previous section from a time serie of absolute returns of the DAX [25] 
during 1973-2002. The most efficient method would be using a maximum like- 
lihood fit. The bad analytical properties of A make this task prohibitive unless 
b is large. As less efficient but tractable method to estimate (3 and r is to use 
X 2 fits to the unconditional pdf given in equ (31). In presence of correlations 
the sample of length T can be still considered as a representative subset of an 
infinite time serie if bT >> 1 holds. A small rate parameter b may cause a 
problem. To see this we divide the data in 5 subsamples of similar statistics. 
Here and in all subsequent estimates we fix r by equ(33) in terms of (5 and 
E'ofv 2 ]. We checked that in each case a free variation of r did not alter the 
results. The estimates of (3 in each sample are given in table 1. They behave 
rather erratic and x 2 probability is rather low indicating a poor descripition 
of the data. This effect is due to the volatility clusters which are by no means 
equally distributed over the subsamples. Taking all data one obtains a good 
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X 2 value and the following value of (3 

(3 = 1.115 ± 0.068 . (42) 

In figure 1 we compare the pdf from equ (31) with the empirical data. The 
agreement is excellent with a x 2 -probability of 60%. The values in the sub- 
samples agree with (42) within two standard deviations. The value (42) cor- 
responds to a tail index of ut ~ 4 which is slightly larger than the the values 
obtained in other financial markets [3,26]. As we discuss lateron this may be 
an effect of not having reached the asymtotic regime in the data. 
We can compare our model with a model discussed by Alfarano [19] where 




0.01 I , , , I , , , I , , , I , , , I 

0.00 0.02 0.04 0.06 0.0S 

v 

Fig. 1. The pdf from equ (31) for equilibrium (solid line) as function of the absolute 
return compared with the DAX from 1973-2002. The dashed line gives the result of 
the linear relation. 

the same herding mechanism is used but the sqare root law replaced by the 
linear relation (11). The fit is worse ( x 2 -P r obability ~ 1% ), but in view of 
possible systematic errors still acceptable. The difference stems mainly from 
the better description of low returns in our model. On a log scale the fit in 
the linear model ( dashed line in figure 1 ) can hardly be distinguished from 
our model, although the tail index tit ~ 6 is 50% larger than in our case. A 
X 2 - or maximum likelihood fit is dominated by the many events with low v. 
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Table 2 

Predicted number of events with normalized absolute return larger than x (column 
1) for the square root law (column 2 and 5) and the linear law (column 3 and 6) for 
two time intervalls. The observed number from DAX is given in column 4 and 7. 

To judge the description of extreme events we compare the predicted events 
with v > x^Je [v 2 ] in both models with the observed number of events. Since 
the largest return in figure 1 correponds to x ~ 6 the fits may not be sensible 
to extreme events with larger x. The values for x > 5 are given in table 2. 
Despite of the large difference in the tail index both models agree with the 
observation reasonably well with a preference for the square root law. Large 
differences occur only outside the observable region (x > 15). We also learn 
that there is no need for an extra mechanism for crash events as proposed in 
[27]. 

Both models can be distinguished if we consider the time dependence. The 
parameter b can be estimated from the common pdf for pairs [vt+t , Vt\ given in 
equ(C5) of the appendix. We perform maximum likelihood fits by adding the 
log likelihoods for all pairs with a time lag t 0i = [1,2,5, 10,20]. Maximizing 
2^j it log f(v t +t i,v t ) with respect to b and j3 we obtain 

P = 1.083 ±0.020 (43) 

b = 0.0075 ± 0.0015 . (44) 



The value of (3 is compatibel with the estimate (42) from the equlibrium pdf. 
The estimate (44) for b implies a large decay time to — 1/& ~ 133[day] in the 
order of 1/2 year. In the linear model only the autocorrelation function can 
be computed. This model cannot describe the data over large range of time 
lags. Adjusting b to describe the autocorrelation for small time lags a value 
b ~ l/37[day]_i is obtained [19]. In figure 2 we compare the autocorrelation 
function given in (39) (solid line) of our model with the values from DAX. The 
agreement is good. The leading asymtotic term (40) (dashed line) shows that 
the data cannot be described by an exponential behaviour. The prediction of 
[19] (dotted line) should be valid only for small time lags. Therefore the data 
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Fig. 2. Correlation function for DAX. Solid line gives the prediction of equ (39), the 
dashed line the asymtotic term (40) and the dotted line the prediction of the linear 
relation. 

for the autocorrelation prefer the square root law. 



Finally we test the scaling prediction, that the conditional pdf f(vt+i\vt) is 
independent of v t if one uses the scaling variable 

x t +i = I — ■ (45) 

y/2/3Eo[v*]+V$ 



The assumption of small b is very well satisfied with (44). In figure 3 we show 
the distribution of x obtained in four intervalls of the normalized absolute 
return 

Vt = i Vt (46) 



at (3 — 0.4,0.7, 1.11,3.0. For a (3 ~ 1.1 the data collapse into a single curve. 
This curve agrees with the predicted scaling function (37) (solid line). From 
this agreement we conclude that the model describes the correlation between 
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Fig. 3. Conditional pdf for the scaling variable Xt+i defined in (45) for various (3 
and vt intervalls. The solid line for (3 = 1.11 gives the prediction 37. The intervalls 
of the normalized return y from equ 46 are denoted by stars (diamonds, triangles, 
squares) for y <0.35 ( 0.35< y < 0.9, 0.9< y < 1.5, y >1.5 ). 

v t and v t +i well, which is not obvious from the autocorrelation due to the 
large fluctuation near t — seen in figure 2. A best fit to the scaling curve 
(37) leads to (3 — 1.1 ± 0.2 which agrees with the previous estimations in (42) 
and (43). It is very remarkable that three very different methods to determine 
(3 lead to the same value. 



6 Conclusions 



The herding model of Kirman has been considered in the limit that funda- 
mentalists are much less affected by herding than noise traders. As a con- 
sequence the market consists of many noise traders and few fundamentalists 
reflecting the conditions in a real market. Volatility cluster occur if the funda- 
mentalists dissapear from the market. This is in contradiction to the results of 
the Lux-Marchesi model, where fundamentalists dominate and volatility clus- 
ters appear if chartists rise beyond 30%. Augmenting the herding model with 
Zhang's square root law between price changes and agent ratios the model 
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turns out to equivalent to the 3/2 model of statistical volatility. This coinci- 
dence allows an interpretation of the parameters in the latter model in terms 
of agent behaviour. Due to the simplicity of our model various observables 
can be calculated analytically. This is important in x 2 or maximum likelihood 
estimations, where numerically accurate and computationally easy accessible 
expressions are required. 

Comparing our model with empirical data we find good agreement with the 
stylized facts derived from the DAX. The tail index of 4 is in harmony with 
universality. Replacing the square root law by a linear relation one gets a 
worse but still acceptable value of x 2 , however the tail index increases to 6. 
This shows that the index depends strongly whether one uses a Pareto law in 
the squared return v 2 (as in our model) or in v. Universality of the tail index 
and the behaviour of extreme events favor the square root law. The autocor- 
relation of the absolute return will decay exponentially at large time lags. The 
very small rate constant b ensures that this asymtotic behaviour is reached 
outside the region of timelags where the data can be trusted. The reasonable 
agreement of the autocorrelation function of our model with the data is much 
better than in the linear model. A quantitave measure of the volatility cluster 
is expressed by a new scaling law for the pdf conditional on the value of the 
previous day, which is again in good agreement with the data. 
The predictions based on a recursive solution of the combined Markov process 
of return and agent ratios will be presented in a future publication. 
Acknowlegments: The author thanks Thomas Lux and Simone Alfarano for 
stimulating discussions and valuable hints. 



A Asymmetric Kirman Model 

In this section we give the changes of the formulae of [18] in the case if one of 
the parameters describing the spontaneous change of oppinion becomes large. 
As pointed out in [18] the macroscopic time scale At during the returns are 
observed does not need to coincide with the microscopic time scale At m over 
which the agents change oppinion. On the latter scale we have the FPE for 
z = n s/( n s + wjr) derived in [18] from the transition probabilities (1) 



For large a\ z will be close to 1. Therefore we introduce the variable u by 



Aw(z,t) _ d 
At m dz 



za 2 — (1 — z)ai + —bz(l — z) w(z, t) 




ai 

u — — 




ai rip 



(A.2) 



b n s + n F 
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which remains finite. Transforming (A.l) into u and neglecting terms of the 
order b/a\ we find 



Aw(u,t) 



a\ ■ D u w(u, t) 



(A.3) 



with the differential operator 



ou 



a 2 , d 
U -J + d~u U 



(A.4) 



The real time scale At (f.e. days) is related to At m by 

At=-±-At m . 
b 



(A.5) 



The limit of large a\ leads to small average values of the number of fundamen- 
talists. By the identification (A.5) the agent ratio can be expressed in terms 
of the time scale 



At r 



n s + rip At 



u 



(A.6) 



Inserting (A.5) into equ. (A.3) leads to 
Aw(u,t) 



At 



= b ■ D u w(u, t) 



(A.7) 



which is equivalent to equ. (7) replacing Aw(u, t)/At by dw(u,t)/dt and a 2 /b 
by (3+1. The FPE (3) ( see [9] ) can be solved by expanding w(u, t) in terms 
of Laguerre polynomials L^{u): 



W(u,t) = J2 Cn(t)L^(u)w (u) 



(A.8) 



n=0 



with the equilibrium pdf w (u) from equ. (7). The functions L^{u) ■ w (u) are 
right eigen functions of D u 



D u ■ L^(m)w (m) = -n ■ L^(m)w (m) . 



(A.9) 



The conditional pdf A(u\u, z) follows from the generating function of the La- 
guerre polynomials (see [29] ). The convolution property (9) is proven by the 
orthogonality of L@ 



E [L^(u)L^(u)} = 6 n , n , 



(A.10) 
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with E [.} the expectation value with respect to wq. 



B Models Outside the Critical Region 



A general class of models can be characterized by the following cost functional 

s(*) = -^E*^ - ^E*. ■ (b.i) 



The minimum coupling go ensures that the system is never in the disordered 
state, uq and /i set the scales for $ and Q. A suitable dynamic (Langevin 
equation, heat bath or Metropolis algorithm) will bring the system into equi- 
librium, which is a Boltzmann distribution. Expectation values are given by 
derivatives of the partition sum Z: 

Z = Uj M*i) exp(-S($) . (B.2) 



For the measure <i/z(<3>) we can use either a spin or a compact model, /i reads 
in these cases 



\ J2$=±a f° r f° r ^e Ising model 
2^ / d$ for for the compact model 



(B.3) 



The expectation value of the return is given by 



v = E 



N v 



o" dlnZ 

N/7 on 



(B.4) 



To evaluate Z we use the Gaus trick for large agent number N: 



<*p(^E*«*i) 



TV 



Nog 



2n(g + g ) 



J dco exp 



2(g + 9o) o"o 



(B.5) 



Inserting equ (B.5) into Z the integration over $j can be carried out with the 
result 
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V 2n{g + ^ ) 
J dco exp ( — TV 



u 2 



(B.6) 



.2(<? + <7o) 

For the spin model we set go = I and C is given by 

C(x) = cosh(x) , (B.7) 

whereas in the compact model we use g = 3 and the function C as 

C(x) = - ■ tanh(x) . (B.8) 

In the large iV-limit Z is given by the maximum of the integrand at uoq. For 
the spin model we get 

= (g + 1) tanh(c<Jo + Hoty • (B-9) 



Replacing ujq by the expectation value v = E[<&i] = uoao/ (g + 1) we find the 
equation of state 

v = a tanh ( (g + 1)— + fj, Q ) (B.10) 



To get the relation with the parameters of the critical model discussed in 
section 3 we expand equ (B.10) around small g and /i 

1 ^ + 1 ^ 3 , 3 - -J-.v - n = o (B.ii) 



3/i V <7 / CT /Xo 

Comparing coefficients with the equation of state (18) we get two equations 
for the three parameters g, <7o and /xo- For example the scale parameter /xo for 
Q can be chosen arbitrarily. Measuring /x in units of Qq we set 

im>= (R -y 3 (B.12) 



and express g and <7o in terms of the previous parameters f2 and as 

5 = o 2 and °o = 7^ 2T ( R13 ) 

3-y 2 a s y{3-y 2 ) 
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Inserting these parameters into equ (B.10) we obtain the equation of state 
valid also outside the critical region 

^•« = 75-^r- tanh (0 (B-i4) 

U (3 - y l )y 



with 



fi 



V + 



V 



6g - 3 



(B.15) 



The parameter y with < y < V3 describes the distance from the critical 
behaviour obtained for y — > 0. The saturation of the return as function of Q 
determines the value of y. The exact form is model dependent. The analogous 
calculation for the compact model leads instead of (B.14) to 

a s _ 45 
Q~o' V= (15-y 2 )y' 

It has the same behaviour for y — * 0, but extrapolates hyperbolically for small 
1/Q to a constant instead of the exponential behaviour obtained from (B.14). 



coth(0 - 



(B.16) 



C Two Point Function of the Probability Density 



From the expansion of the Bessel function in (8) we obtain the following 
representation for A with z = exp(-bt) 



A(u\u, z) = - -— exp 



u 13 ( u zu' 



1 _ z y+i \ 1-z 1-zJ 
00 ( uu'z \ m 1 1 



m=0 



Inserting (C.l) and the noise distribution (28) into equ(34) we obtain after 
integration over u and u' 



f(v v )- N hl-z)^ f ((/? + m + l/2)!; 
m+ * ' Vto) —p\ { } ho m\(J3 + m)\ 



J t+tg 



2rl 



:i + d-^: 



-(/3+m+3/2) 



(C.2) 
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The sum corresponds to the serie expansion of the hypergeometric function 
FM °' X)= E (a-iy. (6-1)! (c + m-1)!^ ■ (C - 3) 



At special arguments F(a, b; c; 0) = 1 and F(a, b; c + 1; 1) = c!(c — a — 6)!/((c — 
a)!(c — 6)!) holds. Using the scaling variable 

C= 2^("«» + < + (1 - z) ^) ( °' 4) 



we get 



'Ar (/5 + l/2)!\ 2 (1-z) 



f(v t +to,Vt ) 



V /?! / (1 + (1 " z)C) W2 
F (p + 3/2, + 3/2; + 1; (1 + ( *_ ^ ) • (C.5) 

With Euler's relation F(a, b; c; x) = (1 - x) c ~ a ~ b F(c - a, c - 6; c; x) equ(C.5) 
can be rewritten as 



F (" 1/2 -" 1/2;/?+1; d + a-»)o ) ' (c ' 6) 

If terms of order 1 — 2 are negligeable and t = 1 ( agrees with the scaling 
variable x 2 +1 from equ (36). Setting in F the argument to 1 leads to 



f{v t+li v t ) = ^±^- {l + x^Y"- 2 . (C.7) 



This equation is equivalent to (35) if we use 

tt i \ dvt+i f{vt+i,v t ) fnQ , 

f(x t+1 \v t ) = - C.8 

dxt+i G (v t ) 



with the equilibrium pdf (31) for v t . 
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D Autocorrelation 



The easiest way to derive E[v% +to ■ vf ] is to perform in equ(34) first the v 
integrations using 

J g{v\u) v« dv = E [v«] {p (D.i) 



with the equilibrium moments 

Ew^. v^tmzim . (D . 2) 



Inserting as in the previous section the expansion (C.I) and performing the 
integrations over u and u' we obtain for the ratio 



R= E[v q t+t0 v q tQ ] 



E [v^ 

e - bt) ^- q . y ( ((3 + m- q /2)\ \ 2 g! ■ 



m=0 



Again the sum is a hypergeometric function leading to 

R = (l- e - u f +1 - q ■ F (j3 + 1 - q/2, /3 + 1- q/2; + 1; e~ bt ) (D.4) 

With the help of Euler's relation we get finally 

R = F(q/2,q/2;{3+l;e- bt ) . (D.5) 

For large bt expansion of F up to linear terms leads to 

R = 1 + 41/TTT) ' e ~ U ' (D ' 6) 



For bt — > and (3 + 1 — q < the expectation value R becomes singular. 
Replacing in equ(D.4) the function F by its value at argument 1 we find the 
power law given in equ (41) 
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